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If the inflaton decays into several components during reheating, and if the corre- 
V^- ■ sponding decay rates are functions of spacetime-dependent quantities, it is possible 

o : 

, to generate entropy perturbations after a stage of single-field inflation. In this paper, 
(N ■ 

' I present a simple toy example that illustrates this possibility. In the example, the 
<D ■ 

ry ! ■ decay rates of the inflaton into "matter" and "radiation" are different functions of 
the total energy density. In particular cases, one can exactly solve the equations 

| of motion both for background and perturbations in the long-wavelength limit, and 

qs, ' show that entropy perturbations do indeed arise. Beyond these specific examples, 

OO : 

(T) , I attempt to identify what are the essential ingredients responsible for the genera- 

(N ; 

t— I ' tion of entropy perturbations after single-field inflation, and to what extent these 

m : 

elements are expected to be present in realistic models. 

I. INTRODUCTION 

& • 
Ctf ■ 

Inflation is part of the modern cosmological standard paradigm because some inflationary 
^ , models can fit the observations . In particular, it is often argued that the "simplest" 
inflationary models [3( yield a homogeneous, isotropic, spatially flat universe slightly per- 
turbed by a nearly scale invariant spectrum of Gaussian, adiabatic density perturbations. 
However, one has to be cautious to generically ascribe to inflation the successes of onlya 
particular set of models. In fact, there are inflationary scenarios that yield anisotropic |^|, 
open or closed [f| universes, strongly tilted spectra 7] and, in multiple-field models, non- 
Gaussian [8( or non-adiabatic [9| perturbations. As it turns out, even single- field inflationary 
models can also produce significant amounts of isocurvature perturbations 1 . 

The simplest inflationary models fall within the class of single-field models. In the latter, 
inflation is driven by a single field, the inflaton, which is assumed to be the only relevant 



1 I use the words "entropy" , "isocurvature" and "non-adiabatic" interchangeably. 
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dynamical component during inflation (in addition to the metric). In particular, pertur- 
bations are imprinted only on the inflaton field. After inflation, these perturbations are 
transferred to the decay products of the inflaton during reheating. Reheating therefore is 
an integral part of any inflationary model. Without reheating, the universe would end up 
dominated by a scalar field in an otherwise empty universe. The reheating mechanism also 
affects the predictions of inflationary models, since observable quantities, like the spectral 
index, depend on the reheating temperature. 

It is believed that, regardless of the details of the reheating stage, the density perturba- 
tions seeded during single-field inflation remain adiabatic. Historically, the first works 
that partially justified this belief assumed that the inflaton decays into a single component. 
If the inflaton completely decays into a single component during reheating, only that com- 
ponent emerges and survives the reheating stage, so no entropy perturbations are possible. 
If the inflaton decays into several components, it has been argued [ll| that perturbations are 
expected to be adiabatic, though I am not aware of a general rigorous proof of adiabaticity 2 . 
In fact, such a proof would be quite difficult to construct, as the couplings of the inflaton 
to matter and the inner workings of the reheating mechanism are hardly ever specified or 
actually known. 

The purpose of this paper is to show that, indeed, perturbations seeded during a stage of 
single-field inflation can be non-adiabatic. In order to do so, it will be enough to construct 
a single counter-example. The essential element of the counter-example is the non-universal 
dependence of the inflaton decay rates on a spacetime-varying quantity, which happens 
to be the total energy density in this particular case. This asymmetry in the decay rates 
finally results in the generation of entropy perturbations from perturbations that are initially 
adiabatic. 



II. A COUNTER-EXAMPLE 



Our goal is to describe the reheating process at the end of single-field inflation and to 
compute the amplitudes of density perturbations at the end of reheating. The basics of 
reheating are reviewed in which focuses on background quantities. In order to follow 



2 See Note added in proof. 
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the evolution of perturbations, additional machinery is needed. Here, I will rely on the 
formalism developed in Il3j], which extends the work of 1141 . This formalism has been 

nnn 

applied in related contexts by several authors 15, 1.16, 1181 ] . 

Suppose that the inflationary stage is driven by the inflaton ip. After the end of inflation 
p oscillates around the minimum of its potential and decays into the particles it couples to. 
In order to describe the decay process it is convenient to resort to a perfect fluid description 
of the different components involved. Let p a and p a generically denote their pressure and 
energy density. These fluids might include the inflaton (a — ip), radiation (a = r) or dark 
matter (a = m). When referring to a generic component, be it the inflaton or its decay 
components, I use greek subindices (a, (3, . . .). When referring to the decay components of 
the inflaton I use latin lower-case ones (a, b, . . .). The inflaton is labeled by ip. 



A. Background 

After the end of a sufficiently long period of scalar-driven inflation the universe is homo- 
geneous, isotropic and spatially flat. In such a universe, the background equations of motion 
of the different fluids during the reheating process are 

pa + 3H(1 + w a )p a = T a p v , (1) 

where a dot denotes a derivative with respect to cosmic time, H is the Hubble constant, 
and w a = p a / p a - The term on the right-hand side of the equation accounts for the decay 
of the inflaton ip into the different matter fields. In particular, — T a is the decay rate of the 
inflaton into the fluid a. Conservation of the total energy-momentum tensor then implies 

^r Q = o. (2) 

a 

It is conventionally assumed that the (scalar) decay rates r^, are constant. Here, we shall 
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In single-field inflation, 



be general and allow spacetime-dependent ones 
spacetime-varying decay rates can arise for instance from a dependence on the inflaton tp, 
or on the temperature of radiation, T a (T), where T ~ p X J^ j^jj. More generally, one could 
consider decay rates that depend on the different energy densities, r Q = T a (pp), which keeps 
the system of Eqs. (JTJ) closed. In the following, I assume for simplicity that the decay rates 
only depend on the total energy density, T a = r a (p). Though this last assumption is not 
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an essential ingredient of the counter-example, it will considerably simplify the analytical 
treatment of the equations. 

If the inflaton potential V around its minimum <p = has the form V oc <p n , on average, 
;he oscillating inflaton behaves as a perfect fluid with equation of state vj v = (n — 2)/(n + 2) 
. For simplicity, I shall assume that all the fluids have the same equation of state w. In 
particular, I assume that the inflaton has the same (constant) equation of state as its decay 
products w = w a . This assumption is also realistic, as we could consider a quartic potential 
for an inflaton that decays into relativistic particles (w = 1/3). 

Summing over all components a in Eq. (JTJ) and using Eq. 0, I obtain the equation of 
motion for the total energy density p + 3H(1 + w)p = 0, which can be readily integrated, 



P = P»H) ■ (3) 

Here, po is the value of the total energy density at an initial time to, which is taken to be 
right at the beginning of reheating. If the decay rates only depend on p, then it is also 
possible to integrate Eq. (JTJ) for the inflaton <p and its decay products a, 

,t / r n 2 



2 

■ (4) 



In the previous equations I have assumed that there are no matter fields at time t , p° a = 0, 
so that p° = po. Because p(t) is known, these equations provide closed exact solutions to 
the equations of motion for any given T a (p). For definiteness, let me at this point assume a 
particular functional dependence of T a . Expand the decay rates in a series in inverse cosmic 
time, and neglect terms of quadratic or higher order. Because the r a depend on the total 
energy density, this translates into 

r Q (p) = r°+ri(£) 1/2 , (5) 

where T° and are two constant coefficients, and po is again the total energy density at 
the initial time t . At any rate, it follows from Eq. (J2J that 

E r ° = E r « = °- ( fi ) 

a a 

The functional form (JSJ) allows me to write the final energy densities in closed form. Plugging 
Eq. (jSJ) into Eqs. (jlj) I find in the limit of large cosmic times 

Pv = Po(f) ^ -exp[r°(t-to)] (7) 
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and 

Pa = p -D. r^ .T(r^ 0) -rJt ,)-^T(r^o + i,-rJto) (y) . (8) 

The (dimensionless) constant D is given by 

D=(-Tlt o y r >.e-r>, (9) 

and T(x,y) is the incomplete Gamma function 3 . Because the inflaton decays, I assume 
that T° and are negative. At late times the evolution of the inflaton energy density is 
determined by the exponential suppression, which rapidly drives p v to zero. Also at late 
times the matter energy density is proportional to t~ 2 , as expected. 

B. Perturbations 

In spatially flat gauge, the perturbed metric has the form 

ds 2 = -(1 + 2<p)dt 2 + 2aB ti dtdx i + a% dx i dx j , (10) 

where <fi and B are scalar metric perturbations. This gauge turns to be convenient because 
in the long- wavelength limit, where one expects spatial divergence of the three- momentum 
and the shear gradient to be negligible |lj|, the perturbation equations do not explicitly 
contain metric perturbations, 

1 S p 

5p a + 3H(1 + w)5p a = --T a p v h ST a p v + T a 5p v . (11) 

1 P 

To arrive at this expression, I have substituted Eq. (37) into Eq. (32) of reference jl3( (note 
that in spatially flat gauge ip = 0). Summing over all components in Eq. ()11|) and using Eq. 
(J2J I obtain dp + 3if (1 + w)5p = 0. Hence, the total density contrast is constant, 

5p = 5_Po (12) 
P Po 

Note that Sp/p being constant does not imply that there are no entropy perturbations. 
Entropy perturbations among components with the same equation of state (the case we 
consider here) do not source changes in Sp/p. Below I will comment on how entropy pertur- 
bations among components with the same equation of state might induce changes in 5p/p 
and thus have observational effects. 



T(x, y) = t x 1 e t dt. I apologize for the idiosyncratic notation. 
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Because the evolution of the background and 5p/p is explicitly known, Eq. (fTTj) can also 
be immediately integrated, 



5p v = 5p 



1 + 



to 



p - -Tp di 



dp r 2^ 



cxp 



r v dt 



to 



and 



Spa 



/ 5p dr a I5p \ ft . 



(13) 



(14) 



't \ Po dp r rv " 2 p ' ' " ' T J \to, 
where I have assumed that there are no matter perturbations initially, 5p® = 0, which by 
the way also implies 5po = 5p° . Substituting Eq. (jSJ) into Eqs. (fT3j) and (fTlj) I get after a 
bit of straightforward algebra 

r£t -2 



i__rj(t_t ) 



■exp [rj(t-to)] 



(15) 



and 



<5p a = 5p - D 



ir° 



5 (rfo - r% - j|) T(i + 1% -r» - ^t(2 + r>, -r>) 



+ (r^o + ^t r>) t(i% -r» 



.(16) 



Note that the Sp a only depend on the T l a through the dimensionless quantity T l a t ~ T % a / H . 
Recall that within the limit of late times, the only approximation that has been made so 
far is the neglect of spatial gradients in the equations of motion (which should be a good 
approximation in the long-wavelength limit) 



Observables 



A set of convenient quantities that characterize the perturbations are the variables ( and 



c 
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23j , which describe curvature perturbations in spatial slices where the corresponding 
energy density is constant, 
5p 1 1 5p 



c 



-H- 



P 



3 1 + w p 



-H 



$Py 
Pv 



H 8 p. 



1 



5 Pa 



and Ca = -#— 

IV p v p a 31 + W p a 



(17) 



The approximations in the last equations apply at late times, when the inflaton has already 
decayed. Entropy perturbations arise whenever the ( a are different for different fluid species. 
A gauge invariant quantity that characterizes such entropy perturbations is 



<S a /3 — 3(Ca — (/?)• 



(18) 
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Entropy perturbations are important because they source changes in (. Assuming that the 
perfect fluids have a constant equation of state parameter, the ( equation of motion is 

C= If^^P^h- (Wa-Wp) -S al3 . (19) 

Therefore, in the absence of entropy perturbations ( is a conserved quantity (recall that we 
only deal here with large scales) . This conservation law has a counterpart for the individual 
( a . It can be also shown that once the inflaton has decayed, ( a remains constant if the 
fluid a is isentropic, i.e. p a = p a (p a )- Therefore, after inflaton decay S a b is constant on large 
scales, regardless of how ( evolves. 

Current experimental results are consistent with a primordial spectrum of adiabatic per- 
turbations |2|, though significant isocurvature components are still allowed j^J. By def- 
inition, perturbations are non-adiabatic if S a p is nonzero for any pair of fluids. Suppose 
momentarily that the decay rates are of the form 

T a = T° a .f(p), (20) 

where T° is a constant and f(p) is any function of the total energy density that is common 
for all species. Because both 5p a and p a in Eqs. ((Hj) and (jU are proportional to T°, this 
constant cancels in the ratio ( a , Eq. (|T7Jl . Consequently, there are no entropy perturbations 
in the matter sector, S a b = 0. In principle, entropy perturbations could arise in the inflaton- 
matter sector, S^a ^ 0. However, if the inflaton decays completely, both its background 
energy density and perturbations vanish, so they do not contribute to the matter budget 
of the universe after reheating; primordial perturbations might then be effectively regarded 
as adiabatic. An alternative interesting — albeit somewhat remote — possibility is an inflaton 
that does not completely decay, say, because its decay rate drops to zero sufficiently fast. 
This frustrated decay might then lead to an inflaton which at late times comes to dominate 
the universe as a dark matter form for instance. Anyway, to summarize: if Eq. (J2*U|) holds 
and the inflaton entirely decays, perturbations are still adiabatic after reheating, S a b = 0. 

If there is more than one component involved in the decay of the inflaton, say, dark matter 
(m) and radiation (r) , and if the decay rates are not constants, there is no reason to expect 
relation (|2()j) to hold. In that case one generically expects non-adiabatic perturbations. The 
source of these non-adiabatic perturbations is the asymmetry in the evolution of the different 
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FIG. 1: A plot of the relative entropy perturbation <S mr /(3£) (z-axis) as a function of T,^ • to [x- 
axis) and ■ to (y-axis). For a large set of decay rates, isocurvature perturbations are significant. 
Note that -t = or T° -t = are ill-defined limits (there is no decay into matter or radiation) . 



decay rates. To illustrate the issue, let me consider the following simple case: 

rj = -r°, r£ = -r^, r°>o, r} = o and r^ = o, r^>o. (21) 

Inserting Eqs. (|2*T|) into Eqs. (JBJ) and (fTBj). and using the definitions (|T7|) and (JTHJ) I find 
that the radiation-matter entropy perturbation is given by 



<? 1 1 S Po 

2 1 + w p 



T(i — r^ 0) r%) t(2 — r^to, r%) 



(22) 



T(-r^ ,r%) T(i-r^ ,r r %) 

Fig. [TJshows the relative entropy perturbation S mr J (3 () as a function of T° t an d r^i . As 
clearly seen in the Fig. ^ the relative entropy perturbation is significant (i.e. of order one) 
for almost any set of decay rates. In particular, in the limit r°t , T^io 1, S mr / (3 () ~ 1/2. 
The sudden drop of the relative entropy perturbation around T^q — is due to the break- 
down of our results for that value (there is no exponential suppression in integrals like (0J)). 

A feature of this concrete example that I expect to survive in models with less assump- 
tions is the linear dependence of the matter density perturbations on the initial density 
perturbation 5p - Hence, if entropy perturbations originate from reheating after single field 
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inflation, the spectral indices of entropy and curvature perturbations should be identical, 
and entropy and curvature modes should be also completely correlated. Though this seems 
plausible and might simply follow from linearity of the equations, as I mentioned at the 
beginning one should be cautious to generalize statements only known to apply in particular 
cases. 

Although our example does indeed contain non-adiabatic perturbations, the Bardeen 
variable ( remains constant on large scales, Eq. (jl2j) . This happens because all components 
have the same equation of state. The last feature though was introduced in the counter- 
example merely to simplify the equations, and it is likely that if the assumption w a = wp is 
dropped while maintaining the non-universal time-dependence of the decay rates, isocurva- 
ture perturbations are still going to be generated during reheating. 

But even if we maintain the assumption of common equation of state during reheating, 
it is possible to generate changes in Q on large scales. Suppose for instance that the inflaton 
decays into radiation (r) and a sufficiently light species of dark matter particles (m). Though 
dark matter might behave as a relativistic component during reheating (w m ~ 1/3), once 
its temperature later drops below its mass, it will start behaving as a pressureless compo- 
nent w m ~ 0. The non-vanishing (and constant) entropy perturbation S mr will then, since 
w m 7^ w r , source changes in ( on large scales and have a significant impact on structure 
formation 251. 



III. CONCLUSION 



In order to show that single-field inflation does not generically produce a spectrum of 
adiabatic perturbations, it is sufficient to present a single counter-example. In this paper, 
I have described a toy reheating process that yields significant amounts of non-adiabatic 
perturbations after single-field inflation. The example contains several ingredients. Some of 
them them are superfluous, and their sole purpose is to simplify the analysis. Others are 
essential, and their presence alone is likely to produce non-adiabatic perturbations during 
reheating. The essential ingredient here is the "non-universal" spacetime dependence of the 
rates at which the inflaton decays into different components. Obviously, such an asymmetric 
dependence can only occur if the inflaton decays into more than one fluid (radiation and dark 
matter for instance). As a caveat, let me mention that I have neglected spatial divergences 
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of the momentum and the shear gradient in the equations of motion for the perturbations. 
These terms are expected to be negligible in the long- wavelength limit, though in order to 
rigorously asses their importance, one should deal with equations of motion containing a 
single variable j^J. 

Is such a non-universal behavior of the decay rates expected to occur in realistic models? 
The universe contains several constituents, like baryons, dark matter, photons, neutrinos 
and quintessence, so there is no reason for the inflaton to decay into a single component. 
In addition, reheating proceeds in different stages. Initially, parametric amplification is 
responsible for an explosive production of particles known as preheating and subsequently 
the inflaton decays by conventional (and essentially different) perturbative processes [12]. 
During this last stage of reheating particles might also acquire thermal masses, which again 
render the inflaton decay rate spacetime-dependent 2lJ . Therefore, it is extremely plausible 
that the inflaton decay rate could depend on spacetime quantities On top of that, 

the decay products of the inflaton might be very different, like dark matter and radiation 
are, so there is no reason to assume that their couplings to the inflaton and, hence, their 
corresponding creation rates are related to each other. Overall, it might very well be that 
the decay rates evolve in a non-universal way. In that case, rather than a spectrum of 
adiabatic perturbations, the signature that points to the single-field origin of the entropy 
perturbations in our example is the complete correlation of entropy and curvature modes, 
and the common spectral index. 

The counter-example presented here does not imply that single-field inflation is in conflict 
with observations. First, significant amounts of isocurvature are still allowed by experiment 



24j . Second, there are models where perturbations still remain adiabatic after reheating. In 



fact, the mechanism I have described can be regarded from a different perspective. Little is 
known about the physics of reheating. Hence, the possibility that isocurvature perturbations 
can be efficiently produced during reheating might, together with observational constraints, 
shed a fair amount of light into our understanding and modeling of the reheating process 
after a stage of single-field inflation. 
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IV. NOTE ADDED IN PROOF 

After this preprint was posted on the archive, S. Weinberg and S. Bashinsky kindly 
pointed out to me that, under quite general assumptions, the adiabaticity of the perturba- 
tions generated during a stage of single- field inflation follows from the results presented in 
27 1 and As suggested by them, the origin of the non-adiabatic perturbations discussed 
here is the discontinuous change of the decay rate along a spatial hypersurface where the 
total energy density is perturbed j^J. Such a change makes the perturbations non-adiabatic 
already at the beginning of reheating. The nature of t he p erturbations generated during a 
stage of single-field inflation is specifically analyzed in 



;he_p( 
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